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We exploit a new theory of duality transformations to construct dual representations of models 
incompatible with traditional duality transformations. Hence we obtain a solution to the long- 
standing problem of non-Abelian dualities that hinges on two key observations: (i) from the point 
of view of dualities, whether the group of symmetries of a model is or is not Abelian is unimportant, 
and (ii) the new theory of dualities that we exploit includes traditional duality transformations, but 
also introduces in a natural way more general transformations. 
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Introduction. — Dualities have been recognized as 
powerful non-perturbative mathematical tools to study 
strongly interacting systems since Kramers and Wannier 
introduced them to determine the exact critical temper- 
ature of the planar Ising model [1]. Traditional dualities 
(TD) as described in Refs. [2-4] are obtained by a sys- 
tematic method based on the Fourier transform (FT), 
suitably generalized to arbitrary groups G. The method 
generates a dual partition function (or lattice Euclidean 
path integral) Z D [K*] from a partition function (PF) 
Z[Ki] with physical couplings Ki, i = 1, ••• , m. The 
dual PF has the remarkable property that its (dual) cou- 
plings K* are large (strong) if the couplings Ki are small 
(weak), and vice versa. This is in part because the dual- 
ity engenders collective (topological) excitations in terms 
of which Z D is expressed. 

Unfortunately, many models of great physical interest 
such as Hciscnberg, non-Abelian gauge and more recent 
models based on Hopf algebras are outside the scope of 
TD transformations. The reason is technical, not physi- 
cal: the group-theoretic FT has different algebraic prop- 
erties depending on G being Abelian or not, and the 
TD transformation takes advantage of essential simpli- 
fications present only in the Abelian case. In essence, a 
TD transformation introduces, via an FT, dual elemen- 
tary degrees of freedom (EDFs). For Abelian FTs, the 
dual EDFs arc still locally coupled and result in physical 
dual PFs. Non-Abelian FTs result in non-local interac- 
tions and/or constraints and complex Boltzman weights, 
as historically illustrated by attempts to construct dual 
representations of non-Abelian gauge theories [5]. Thus, 
in order to obtain TDs, it is necessary that the model and 
associate groups satisfy restrictive properties enabling 
the existence of physical dual models. 

Conventionally it is thought that the group G needed 
for TD transformations is determined by the model's 
group of symmetries Q [2, 6] (see especially section 7, 
point (3) of Ref. [6]). Here we argue that G is not deter- 
mined by, and in general is unrelated to, Q. Rather, G 



is associated with and constrained by the model's local 
or quasi-local interactions. We call a model S-Abelian, or 
S -non- Abelian, according to whether the group of sym- 
metries G is Abelian or not. Many models are S-non- 
Abelian, but have a TD transformation with an associ- 
ated Abelian G. It is tempting to call a candidate duality 
transformation D- Abelian or D-non- Abelian according to 
whether G is Abelian or not. However, the underlying 
group may not be apparent and may involve more gen- 
eral structures. Instead, we focus on the presence or 
absence of non-trivial constraints on the states of the 
models. That is, we say that a transformation connect- 
ing two locally defined PFs has D -non- Abelian features 
if the transformation introduces or removes non-trivial 
local constraints. From this perspective, it is impossible 
to have a D-non- Abelian self-duality. 

The non-Abelian duality problem is the problem of ex- 
tending the scope of TDs without sacrificing their physi- 
cal content to cases where there are no relevant Abelian 
groups G for the interactions of a model. Our main con- 
tribution is to introduce a generalization of TD transfor- 
mations, bond- algebraic duality transformations, that ad- 
dresses the problem of non-Abelian dualities by exploit- 
ing the recently developed theory of bond algebras [7, 8] 
and their homomorphisms. These transformations [2] 
handle on equal footing models with arbitrary G, Abelian 
or not, and even more general models, where there is 
no obvious group structure constraining the transfor- 
mations. Unlike a strictly D-Abelian duality, a bond- 
algebraic duality can have both D-Abclian and D-non- 
Abelian features. To illustrate our ideas, we give a du- 
ality for a model outside the scope of TDs, namely a 
rigid-rotator model with group G = SU(2). According 
to our terminology, this duality is D-non-Abelian and 
impossible to obtain by a TD. 

Lattice Models. — For simplicity, consider models with 
identical, classical EDFs with configuration space M 
at sites r of a lattice A. A full configuration of the 
model consists of an assignment s r € M for each site 
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r. If the model has only pair- wise symmetric interac- 
tions, then the total energy E{s r } of a configuration 
{s r } is a sum of (oriented) two-body interaction energies 
e(s r ,s r i) = e(s r i,s r ). This minimal description suffices 
to specify physical quantities such as a PF. However, it 
often happens that M admits useful additional mathe- 
matical structures. In the context of TDs, this includes 
groups acting on the EDFs. More generally, we can con- 
sider configuration spaces that are endowed with two op- 
erations m, ml i-> m • m! (multiplication) and to t— > S(M) 
(involution) such that (a) multiplication is associative, 
(b) S is involutive (<S 2 is the identity map) and order- 
reversing (S'(m-m') = S(m')-S(m)), and c) the pair-wise 
interactions between EDFs can be expressed in the form 

e(s r , s r ') = v(s r ■ S(s r ')*\, (1) 

for some real- valued function v. Conditions (a) and 
(b) turn M into a semigroup with involution. We call 
models satisfying these conditions m-models (short for 
multiplication- models) . It is possible to accommodate in- 
teractions involving more than two EDFs, provided the 
EDFs in an interaction are ordered and oriented. For 
example, let s ri , s r2 , s r3 , s ri occupy the corners of an 
elementary plaquette on the lattice, ordered along the 
boundary of the plaquette. Then 

^(^t*i j s r2 , s r3 , s r4 ) v (s ri • 5*(s r2 ) • s r3 • S(s r4 )) (2) 

describes a form of m-interaction relevant to physical ap- 
plications that we discuss in the next section. 

Wilson's lattice approach to quantum field theory [10] 
popularized the study of m-models defined in terms of 
EDFs taking values on a group G = M, with interac- 
tions of the form of Eq. (1) or its generalizations. These 
G-models are important examples of m-models where 
the multiplication in M is group multiplication and S 
is group inversion, S(g) = g . TD transformations 
are applicable only to G-models with G an Abelian group 
[3] . A reason for introducing the more general notion of 
m-model is that we want to accommodate a larger set 
of theories, such as those based on general Hopf algebras 
[11] that are becoming increasingly more important in 
topological quantum matter, and the theory of quantum 
computation and error correction. 

A model's symmetry group Q is completely determined 
by its interactions. But semigroups with involution M as- 
sociated with the model and constrained to satisfy iden- 
tities such as those of Eqs. (1) or (2) are in general not 
unique and may be completely unrelated to Q . For exam- 
ple, consider the non-Abelian group Sn of permutations 
on N > 3 letters, and use it as the configuration space 
AI = Sn for the EDFs of the Potts model. Then we can 
write the interaction energy as 

ePotts(Sr, S r ') = S e (s r ■ s^, 1 ), (3) 

where S e (g) — <5 e , s is the Kronecker delta on Sn- The 
Potts model is non-Abelian from the point of view of its 



symmetries, but it supports D-Abclian dualities. The 
reason is that we can map the elements of Sn to the ele- 
ments of Zn\ (the Abelian group of integers modulo AH), 
and rewrite the interaction energy in the equivalent form 
£potts(sr, s r >) = <5o(s r — s r >). Rewriting the model in this 
way does not change the fact that its symmetries are non- 
Abelian, yet it permits the use of a TD to determine its 
critical coupling. Some early explorations of non-Abelian 
dualities [12-14] exploited this procedure extensively to 
map models defined on certain non-Abelian groups to 
Abelian ones. In particular, it was noted that models 
defined on solvable groups are specially amenable to this 
procedure [13], since solvable groups can be mapped to 
Abelian groups in a natural way. 

Beyond traditional dualities. — The recently developed 
theory of bond-algebra homomorphisms [2, 8] includes 
and generalizes the theory of TD transformations. To 
apply this theory, we start with a physical model de- 
fined by its EDFs and local interactions that capture the 
main features of the physical phenomena under study. 
We then identify the model's bonds, which are the local 
or quasi-local interaction operators occuring in the inter- 
actions. The multiplicatively closed algebra generated by 
the bonds is called the bond algebra. A key observation 
is that the structure of the bond algebra and its generat- 
ing bonds contain essential information about the model. 
In particular, mappings between bond algebras that pre- 
serve locality in, and all the algebraic relations among 
the bonds, can demonstrate close relationships between 
seemingly unrelated models, including models with EDFs 
of differing exchange statistics. Although such bond- 
algebra mappings are by definition local in the bonds, 
they are typically non-local in the EDFs. That is, the 
model's EDFs in the domain can be naturally related 
in the range to highly non-local degrees of freedom in- 
volving many EDFs [2, 8]. These collective modes can 
be considered to be alternative EDFs relative to which 
interactions take different, but still local, forms. In the 
following, we call mappings of bond algebras that pre- 
serve locality and algebraic relationships bond- algebraic 
duality transformations. This is motivated by the obser- 
vation made in Rcfs. [2, 8] that they can be used as the 
foundation for a unified theory of classical and quantum 
dualities. Here we show that bond-algebraic dualities go 
beyond TDs and generate new transformations that are 
not related to the group-theoretic FT. 

A bond-algebraic duality [2] for a classical model can 
be obtained by expressing the PF Z in terms of opera- 
tors that can be related to a bond algebra. A popular 
way to do this (for an alternative, see [15]) begins by 
identifying operators Tq, ■ ■ ■ ,T S , called transfer matrices 
(TMs), acting on a Hilbcrt space H, and a preferred basis 
<f> — {\4>i)} of W. The operators must satisfy 

2 = E le~ E{ST} ] = Tr*[(T. ■ ■ ■ TiT ) N ], (4) 
{ St .} 
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where N is determined by the length of the lattice in a 
chosen direction. The role of the basis is so enable us to 
make the equality explicit by appropriately inserting res- 
olutions of the identity |<fo)(0i| between the operators 
in the trace, expanding the trace in terms of the resulting 
summands and associating the states s r with sequences 
of basis indices. For this to work and the expanded trace 
to match the desired PF, we need the right combination 
of TMs and a preferred basis. 

The locality of the classical model's interactions is usu- 
ally reflected in this construction. Thus, the Hilbcrt 
space H is defined by quantum EDFs on a lattice such 
that the TMs factor into a product of quasi-local oper- 
ators, T a = Yir^a,r ( a = 0, l,--- , s), with T a lattice 
index that may stand for a site, a link, or a plaquette. 
As a result, it is natural to define the bond algebra of Z 
as the algebra generated by the bonds {t a ,r} [7]. 

To obtain a duality, one can algebraically represent the 
bonds, and therefore the TMs, on an alternative space, 
and determine a preferred basis so that the expansion of 
the trace can be recognized as a physical PF for a differ- 
ent model. Suppose we have such a bond-algebraic dual- 
ity t Qj r ^ i„p with image bonds t% r on different quan- 
tum EDFs that are also local and have the same algebraic 
relationships. This induces a bond-algebra isomorphism 
between the algebras generated by the two sets of bonds. 
We can define dual TMs T£ = A~ 1/N ]J r t£ r , with A a 
analytic functions of the parameters of the model, and 
compute a dual PF as 

Z D = Tr4(T s D -..T 1 D T D ) N }. (5) 

relative to a basis ip = to b° specified. A non- 

trivial property of typical bond algebra isomorphisms is 
that they are induced by unitary transformations [2]. In 
particular, if r = Uit a jU^ with unitary, then 

s 

Z = AZ D 8M& A=Y\A a . (6) 

Q=0 

It follows that Z and Z D represent two, in general dif- 
ferent, systems that have nonetheless the same thermo- 
dynamics. 




FIG. 1. Lattice connectivity of the classical D = 3 XM model. 

The final form of Z D in terms of its EDFs depends 
critically on the choice of basis ip in Eq. (5). As an 



extreme example, if {\ipj)} — {Ud\(j)i)} } then Eq. (6) is 
reduced to a trivial identity with A = 1. The choice 
of basis also determines whether a bond-algebraic dual- 
ity is D-non-Abelian or D-Abelian, that is, whether or 
not it introduces local constraints when the trace is ex- 
panded. Local constraints appear if the combination of 
TMs between resolutions of the identity have entries that 
are zero with respect to the basis. Thus, given a bond- 
algebraic duality, it is natural to seek a basis where the 
relevant TMs are full, so that the duality is D-Abelian. 
In general, the entries of the matrices also need to be 
positive and expressible as products of local Boltzmann 
weights. Although such bases are known to exist for a 
large class of duality problems including TDs, we do not 
have general strategies for finding them. 

We illustrate these ideas with a D-Abelian and a D- 
non-Abelian duality for the Xu-Moore (XM) model of 
p+ip superconducting arrays [16, 17]. The model's D = 3 
dimensional classical PF is given by (see Fig. 1) 

Zxm = e E s E T .[A"i(T (T ., i)C r (l ,, I+1) +A' p n CT ( T ,. i) ] i ^ 

where cr^ ri ^ = ±1 are classical Ising variables placed 
at the sites (r,i) (i an integer) of a cubic lattice, and 

D^r.i) = V(r,i)V(r+e 2 ,i)0~{r+e 1 ,i)V(r+e 1 +e 2 ,i)- The XM 

model is a G- model with G = Z2. The TD transfor- 
mation maps the model to itself with a characteristic 
interchange of strong and weak coupling constants [16]. 
To recast it as a D-Abclian bond-algebraic duality, we 
construct plane-to-plane TMs 

T 1= n^, To=n eW . (8) 

r r 

with cr£' 2 Pauli matrices acting on quantum spins at sites 
r of a [d — 2) square lattice, h = — lntanh(_ftT/)/2, and 
Ua% = cr^ +e2 a^ +ei a^. +ei+e2 (see Fig. 2). To recover 
Zxm, the trace Tr^TiTb)^] is computed with respect to 
the basis <f> that diagonalizes the a%. 

The TMs can be expressed as products of t\_ r = 
cosh(/i)+sinh(/i)cr*,t 0iT , = cosh(7^ p ) + sinh(_ftT p )D(7 2 . We 
therefore let the bonds be {a^,Da^.}. They satisfy a 
bond-algebraic duality induced by 

al^Ual, Ua z r k> <r* +ei+ea , (9) 
and illustrated in Fig. 2. The dual TMs 

T? = I] e hD< , T a D = [] e K *< , (10) 

r r 

are related to T\, Tq by a unitary mapping. If we expand 
z xm = Tr 4>[{ T F T o') N } wi th respect to fj>, we find that 
^xm contains local constraints, so that the mapping of 
Eq. (9) is D-non-Abelian relative to (j>. It is, however, 
£>-Abelian in the basis ifj that diagonalizes the <j%, with 
respect to which we recover the traditional self-duality of 
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FIG. 2. The quantum XM model (shown on top) is self-dual 
as indicated by the arrow on the left, and it is dual to the 
planar orbital compass model, as indicated by the arrow on 
the right. Direct and dual lattices are indicated with solid 
and dashed lines, respectively. 



the XM model [8, 16]. In the bond-algebraic approach to 
dualities, the role of the FT is encoded in the change of 
basis (j) ' — ^ tp realized by a direct product of Hadamard 
operators H satisfying Ha z H = a x . 

The bond algebra of the XM model has another local 
representation [2, 8, 18], 



r+e 2 ' 



□erf 



The corresponding dual TMs 



(11) 



(12) 



-"XM 



yield an alternative dual partition function Zy 
Tr ^[(ffff ) N ]. With <t> the basis that diagonalizes the 
<j% , we obtain a PF with local, four-spin constraints. Rel- 
ative to this basis the duality of Eq. (11) is D-non- 
Abelian. It is an open problem whether there is a choice 
of basis for which Z® M is an unconstrained canonical en- 
semble making the duality D-Abelian. An alternative 
may be to remove these constraints by reinterpreting 
them as gauge symmetries. 

It is important to recall at this point that a TD maps 
a G-model on a lattice A to an essentially unique dual 
model supported on the dual lattice A* [3], and the XM 
model is self-dual under such TDs. In contrast, the bond- 
algebraic duality of Eq. (11) results in a model with a 
Hamiltonian that differs from that of the XM model. We 
conclude that this bond-algebraic duality is not a TD. 

Non-Abelian dualities. — Next, we show that bond- 
algebraic dualities exist for G-modcls with non-Abelian 
G and no TDs. For example, consider the Euclidean 



lattice version [10] of the SU(2) principal chiral field 
[19]. This model involves an S'/7(2)-valued field u(x) = 



G SU(2) with action 



Spa = 77T — / dtdx tr(d u* .d u — d\u* .diu) 
2A J 



(13) 



The lower dot denotes matrix multiplication, u* (x) is the 
Hermitian-conjugate field, and tr is the 2 x 2-matrix trace. 
Since u*(x)u{x) = 1, the lattice Euclidean path integral 
reduces to 



PCh 



Re{tr(ti; +ei « T .)+tr(ti; +e2 « r )} 



{«-r} 



(14) 



on the square lattice with SU{2) as the EDFs' config- 
uration space. Note that if we replace SU(2) by U(l) 
we obtain the AF-model, for which there is a D-Abelian 
duality to the solid-on-solid model [2]. 

To express Zpa in terms of row-to-row transfer oper- 
ators, we use covariant pairs of standard representations 
of SU{2) and the continuous functions C$(SU(2)) on 
SU(2), both acting on wavefunctions on SU(2). A gen- 
erating set for C (SU(2)) is given by (U)» u {fx, v = 1, 2), 
where ([/^(it) = u tl v . Thus U is a matrix-valued func- 
tion. The standard representation of SU (2) has infinites- 
imal generators J = {J Xl J y ,J z ) for multiplication on the 
right. If we write u = e - |6,n tJ / 2 j Q a finite angle, and n 
a unit vector, then e' 9n ' J \v) = \v.u) for the formal basis 
of wavefunctions \v). The row-to-row transfer operators 
are given by 



T =n 

i 

Ti = n / du 



e ^Re{tr(C/*.C/ I+1 )} 



e fcRe{tr(«)} e iflft.Ji u = e 



(15) 



J ia — 

0. The 
Ji and 



for a parameter h dependent on Ao. The products are 
over the EDFs in a row. To recover Eq. (14), the trace 
Tr^[(TiTo) w ] is expanded with respect to the basis \v). 

To define a bond-algebraic duality, we use the genera- 
tors ji of left multiplication, which satisfy e = 
\u.ui). These generators can be related to actions 
defined by J and U by the identity [20] 
Eb=x., v ., z tr ( u ^ a U l <J b )J lb /2, such that [j^Jj] = 
bond algebra generated by the local bonds 
U* .Ui+i can be transformed to local bonds according to 

Ji h> -ji + Ji-i, U*.U l+ i i ^ Ui. (17) 

Proving that the mapping is induced by a unitary op- 
erator requires adding boundary terms to complete the 
algebra, checking that the images of the EDFs' operators 
are generated by corresponding covariant pairs of repre- 
sentations and applying the Stone- von Neumann-Mackcy 
theorem [21] (see the Supplemental Material). It follows 
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that 



T Q D = H e^Mc/*)} t (is) 

i 

T ? = II / du e x ^< u ^e m <~ M -^ (19) 

are unitarily equivalent to the corresponding TMs To and 
T\. Note that the dual variables Ji, Ui 



J' i ji ~\~ Ji- 



Ui = --- .U* +2 .U* +V U* (20) 



that are the unitary images of the EDF operators under 
the duality are, as expected on general grounds, non- 
local collective modes. The string defining Ui extends 
to the boundary of the system, and its specific form is 
determined by the chosen boundary conditions. 

To obtain a dual PF, we expand the trace with re- 
spect to the basis \vi) for each i. The PF Zp Ch = 
Tr^,[(TfT^) N ] is then given by (r = ie ± + je 2 ) 



f e E« 5^ Re { tr «,j-»i,i+0+%,2j)} Y[ 6(1, Oi 

{"*•} i+j=even 



i), 



where Ouij = u* jUij + iUi+i j+iu* +1 j (see the Supple- 
mental Material). As for Z® M , we obtain a PF with local 
constraints on a checkerboard. We do not know whether 
there is a choice of basis that removes such constraints. 

We have discussed dualities for classical models, which 
also apply to Euclidean path-integral representations 
of quantum problems. This is the context in which 
the problem of non-Abelian dualities is typically stated. 
However, as explained in detail in Ref. [2], bond- 
algebraic dualities provide a unified approach to classical 
and quantum dualities, so we can use essentially the same 
techniques to obtain dualities for any quantum mechani- 
cal model. The bond algebra of a quantum Hamiltonian 
H = J2r ^ r * s t ne algebra generated by the local or quasi- 
local bonds hr, and a bond-algebraic quantum duality is 
given by a mapping hr to an algebraically equiva- 

lent dual set of local or quasi-local bonds. As before, one 
can typically show that the isomorphism is induced by a 
unitary transformation, in which case H D = hp is 
unitarily equivalent to H. Take, for example, the d = 1, 
infinite chain, SU (2) equivalent of the Z 2 transverse-field 
Ising Hamiltonian 



A 



which is not self-dual, but has a duality to 
tf&h = + ^-i) 2 + £(f(E?) + tr(Di))], (22) 



so the distinction between D-Abelian and D-non-Abelian 
becomes irrelevant. 

Contributions to this work by NIST, an agency of the 
US government, are not subject to copyright laws. 
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as follows from Eq. (17) (see the Supplemental Material). 
Quantum dualities are remarkably simpler than classical 
dualities. They do not depend on a choice of basis, and 
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SUPPLEMENTAL MATERIAL TO THE 

NON-ABELIAN DUALITY PROBLEM 

We present mathematical details and clarify technical 
issues of results reported in the accompanying paper The 
non-Abelian Duality Problem. 



THE ISING MODEL REVISITED 

In this section we explain how to choose the basis 
for expanding the traces when determining partition 
functions from products of transfer matrices. We illus- 
trate the main concepts by example and use the two- 
dimensional Ising model [1] to show that bond-algebraic 
dualities may display both D-Abelian and D-non-Abelian 
features. The partition function of the Ising model is 
given by (r = ie x + je 2 ) 

Z\[K 1 ,K 2 ] = V] exp[y^(.K'icri i jcr i+1 j + K 2 aija iij+ i)}, 



{<M 



i,3 



(23) 

and can be expressed as Z\[Ki,K 2 ] = Tr^[(TiTo) N ] in 
terms of the transfer matrices 



r =n ^ c 



T^Hie^ + e-^af), (24) 



provided the trace is expanded in the basis <p = {|0fe)} 
that diagonalizes the Pauli matrices of. 
The mapping 



a*a z l+1 h-> of , af^af_ t af 



(25) 



illustrated in Fig. 3 defines an isomorphism of bond alge- 
bras that is induced by a unitary mapping. Thus 7b, Ti 
are dual and unitarily equivalent to 

T D = [] e K ^, T? = J] {e K '+e- K 'afaf +1 ). (26) 



X 

a i-l af 1( T 




FIG. 3. Duality isomorphism of bond algebras associated with 
the transfer matrices of the Ising model. 

To compute a partition function from the dual trans- 
fer matrices via the expression Zf = Tr^[(T-j D 7 , c P) JV ], we 
need to specify a basis ip = {\ipk)}- The expansion of 
the trace obtained by inserting resolutions of the iden- 
tity with respect to this basis must be recognizable as 



the partition function of a local system. In particular, 
the coefficients of the expansion must be non-negative, 
so that they can be written as Boltzman weights, and 
they must be products of local terms consistent with the 
expansion. For example, set ip = <f>, the basis of the pre- 
vious paragraph. Then, as will become clear below, it is 
convenient to split Tf = Tf D Tf e , where 



T 



D 



e"^« +1 ), 



K-2 



-K2 „x „x \ 



(27) 



(28) 



with i Q = 2i+l, i e = 2i, i 6 Z. We can label the members 
of the basis \(f>k) in terms of strings a of Ising variables 
±1 at sites i so that af \a) = a- t \a) . With these labels, 
the basis members are written as \a). We now compute 
Z F = Tr 4>[( T Fo T Fe T o') N } b y expanding the trace as 

T^[(TKT„) W ] = (29) 
]T {a l \TE ) \a 2 ){a 2 \T^\a 3 )-.. 

{<Ti},— ,{(T2Jv} 

' • ' (o'2Ar-l|7 , f^|(T2Ar}(o'2Af|7 1 i^7o D | Cr l)' 

where = 1, • ■ • , 2N, describes the state of row j. 

Note that Tjf is diagonal in the chosen basis. Further, 



( a ]\ T Fo\ a j+i) 



(30) 



n 



K 2 



and a similar factorization holds for (aj\T-[^\aj + i). The 
splitting Tj° = T^T^ was introduced to ensure this fac- 
torization. We can evaluate Eq. (29) by applying shifted 
forms of the identity 



—K2 x x 



e 2 



(a^ai+a^ 02) 



a\a 2 \a\a 2 ) = 
S(a' 1 ai,a' 2 a 2 ) 



(31) 



It follows that 



{ay} \i+j=even 

x cxp (-^-Vi,i<ri,j+i + Kia h2j+ i] 



(32) 



h3 



The interactions in Zf are illustrated in Fig. 4. 

The last factor in Eq. (32) for Zf can be identified 
as a Boltzmann weight for a physical system with local 
interactions, one of the requirements for a good choice 
of basis to expand the trace in. However, the expression 
for the partition function in Eq. (32) also introduces lo- 
cal (delta function) constraints to account for the fact 
that the dual Boltzmann weights vanish for some config- 
urations. It is preferable to find a basis ijj where all the 
Boltzmann weights are strictly positive so that there are 
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x x x 



from Eq. (25). Hence 



"i — i — i — r 

FIG. 4. Interactions and constraints in the dual partition 
function Zf . The crosses highlight the sites where the classi- 
cal Ising variables couple to a inhomogeneous external field of 
magnitude K\. The heavy vertical lines indicates a nearest- 
neighbor Ising interaction of magnitude Kijl. The staggered 
distribution of plaquettes with round corners indicates the 
distribution of four-spin delta constraints. 



m,n&m-\- X,n^Tn-\-2,n 



(35) 



In the dual model Zf, the string correlator (/i TO ' >n ' ix m ,n) 
is (in the limit of infinite separation) the square of the 
order parameter. Thus, for example, if K < K c , Z\ is 
in its ferromagnetic phase, corresponding by duality to a 
phase of Z® dominated by strong correlations of string 
collective modes. 



DUALITY OF THE SU(2) PRINCIPAL CHIRAL 
FIELD: HAMILTONIAN FORMULATION 



This section discusses a duality for the finite system 



N N-l 

Hpcu = 2 E J ™ + 2 E Retr(£C +1 .(7 m ). (36) 

m— 1 m— 1 



no constraints. For the Ising model, one can find such a 
basis by inspection. Let ip be the basis that diagonalizcs 
the Pauli matrices af . Then one can check that 



The Hamiltonian iJpch can be obtained as the time- 
continuum limit [3, 4] of the partition function of Eq. 
(14) of the accompanying paper. 



Z\[K U K 2 



Tr,[mT ) w ] 



AZ\[Kl,K* 2 ] (33) 



Algebra of a Single Quantum Rigid Rotator 



with sinh(2i^ 1 *)sinh(2A' 2 ) = 1 = sinh(2/q) sinh(2ii:i). 
The proportionality factor A is an analytic function of 
the couplings and size of the system [2] . We then recover 
the Kramers- Wannier self-duality of the Ising model. 

The duality of the Ising model expressed by Eq. (32) 
is not a self-duality. The Kramers- Wannier self-duality 
as derived above is the result of combining the bond- 
algebraic mapping of Eq. (25) with a suitable choice 
of basis tp. The dual partition function according to 
Eq. (32) has restructured the interactions drastically, 
but has left the couplings K\ , K 2 essentially unchanged. 
Nevertheless, such dualities reveal key properties of tra- 
ditional dualities. For example, consider the two-point 
correlator (er m ',„' cr m n ). In the limit in which (m',n') 
is infinitely far from (ro,n), this correlator defines the 
square of the order parameter. We can compute the cor- 
relator in the dual model of Eq. (32) as 



Tr^r^-"') a z m , T("'~") a z m T n ] 

Tr<j>[{T D ){ N-n') ^ {TD) (n'-n) ^ {T D )n] 



(34) 



N] 



(fJ>m' ,n' f^m,n) •> 



where T = TiT , T D = T®T® (sec Eq. (27)), and a z H> 
m fT m+i (J m+2 '" > the dual Pauli spin operator 



The kincmatical algebra of a rigid rotator [5] is defined 
by the relations among the canonical variables J a , a = 
x,V,z, = 1,2, 



4 

[Ja, Jb] 



J a i 

^abcJ c 

l 



[J,U] = -U.a, 
U*.U = U.U* = 1, 



(37) 
(38) 

(39) 
(40) 



introduced in the accompanying paper. Here a denotes a 
standard Pauli matrix. The low dot denotes matrix mul- 
tiplication to distinguish it from tensor multiplication, 
and a centered dot denotes the standard Euclidean inner 
product. For example, J • J = J% + J? + J%, and 



[J,U] 



l -U.a O [J a ,U\] = \ 



(41) 



Eqs. (37) and (40) imply [J, U*] = -a.U*/2. 

The algebra above affords a set of position-like oper- 
ators U, U* and conjugate momenta J a that suffice to 
specify completely the kinematics of quantum tops. It is 
useful however to introduce three additional operators 



Ja 



^J2tr(U*.a a .U.a b )J b , 



x,y,z 



(42) 
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or just j = tr(U* .a.U.(a ■ J))/2 for short, having some 
very useful properties: 



J a 3a i 
[jaijb] — i^abcjc: 

\j,U] = ±a.U, 

[ja,Jb] = 0, 

3-3 = J- J- 



(43) 
(44) 

(45) 

(46) 
(47) 



Direct proofs of these relations, based on definition Eq. 

(42) and relations (37), (38), (39), (40), can be found in 
Sect, of this Supplemental Material. Notice that Eqs. 

(43) and (45) imply that \j,U*] = -U*.a/2. 



Bond-algebraic Duality Transformation 

This section describes the construction of a dual rep- 
resentation of the Hamiltonian Hpch of Eq. (36). The 
starting point is the selection of a suitable set of bonds 
as generators of the bond algebra of interactions. One 
convenient choice is 



J n 



m = 1, 



,N 



(48) 



u; n+1 .u m , u; n .u m+ll m = 1, • • • , n - 1, (49) 

We call the algebra they generate *4pch- Notice that 
Hpch G -4pch, but the bond algebra does not include 
the position-like operators U m ,U* nl m = ,N. It 

will be useful later to change this by adding a boundary 
term 



U N , 



u* N , 



(50) 



to the list of generators of Apcb- The resulting ex- 
tended algebra, still denoted by Apcu, does include the 
1, ■ • • ,N, since 



U m ,U^, m 



U m = U N .(U* N .U N ^). ■ ■ ■ .(U m+1 .U m ), (51) 
U* m = {U* m .U m+1 ). ■ ■ ■ .(UZ_ 1 .U N ).UZ . (52) 

The extended algebra Apa is simply a direct product 
of N copies of the algebra generated by a single rigid rota- 
tor J, U, U* . However, what is required is an understand- 
ing of the structure of *4pch from the point of view of the 
local interaction terms in iJpch • The relations (other than 
commutation) between the bond generators of Eqs. (48), 
(49), and (50) are U* N .U N = 1, {U^ +l .U m ).{U* m .U m+1 ) = 

1, [J m ,a, Jn,b] = itabcJm,cbm,n for 771 = 2, • • • ,N — 1, 

[J m ,U* n .U m+ i] = —-(T.U^.Um+i, (53) 

[Jm, U* m -M = \u* m _ x .U m .a , (54) 

and at the boundaries, [Ji, U^.U^ = —\a.U*.U2, 
[Jn,Uh_ v U n ] = ^Uff^.UN.a, and [J N ,U^] = 




FIG. 5. Duality automorphism for the quantum chain of rigid 
rotators, shown for three sites (N = 3). 



-2 a - U N- 



Relations that follow by Hermitian conjuga- 
tion from those listed have been omitted. 

The goal is to construct a mapping that preserves these 
algebraic relations and locality. For instance (see Fig. 5), 

U* N ^U N , CAh^i, (55) 

Jl >->• -ji, J m >-> -jm + Jm-1, (56) 

for m = 2, • • • , JV. It is not necessary to specify the action 
of this mapping on the j m , since the j m are functions of 
J mi U m ,U^ (secEq. (42)). 

As noted in the accompanying paper, to verify that 
the bond-algebra mapping defined above is induced by 
a unitary map, we can invoke the Stone-von Ncumann- 
Mackey theorem [6]. In order to do so, we need to verify 
that the operators of the elementary degrees of freedom 
are transformed into operators of a covariant pair of rep- 
resentations as required by the theorem. We can express 
the images of the operators in terms of the bonds (in- 
cluding the boundary terms) directly. A benefit of doing 
so is that these images define collective modes of interest. 

The dual momenta J m are by definition the image 
Jm !- > J-nn and are obtained directly from Eq. (56), 



Ji 



J\ It 



Jrn — l jr 



(57) 



for m = 2, • • • , N. To compute the dual position-like 
operators it is necessary to exploit the decompositions 
of Eqs. (51) and (52). These decompositions combined 
with Eq. (55) yield 



U m ^ U 



N ■ 



■ U* = u n 



(58) 



and U* % i-> U m . ■ ■ ■ .Un = U m - It can be checked that the 
dual variables Jm.Um^U^ commute on different sites, 
and satisfy the relations of Eqs. (37), (38), (39), and 
(40), as required for a covariant pair of representations. 

Similarly to the dual variables, the dual Hamiltonian 
is computed as £Zpch 



Hp Ch . Hence 



JV-l 



A 



+i- JmY + -Retr(C/ m ) .(59) 



m— 1 



To gain insight into the physical meaning of Eq. (59) 
it is useful to discuss the global symmetries of -Hpch and 
their dual representation. On one hand, the interaction 
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terms Retr(U* n+1 .U m ), m — 1, ••• ,N — 1, arc invari- 
ant under right and left multiplication, U m — > U m .v 
and U m — > v.U m . It follows that iipch has a global 
SU(2) x SU(2) symmetry, with infinitesimal generators 

^ = Em=i^ and = El=iim tnat commute with 
i/pch- On the other hand, the dual Hamiltonian i/pch 
contains the terms Retr(C/ m ), m = 1, • • • , N — 1, invari- 
ant only under the adjoint (anti)action, tr(v* .U m .v) = 
tr(f/ m ). It may seem that a symmetry has been lost. 

The duality maps the symmetry generators j, J to dual 
symmetry generators 

N N 

J = --/m = -jl + ^ (~3 m + J in-l), (60) 

m— 1 m— 2 

AT N I 

J = = E 2 tr (^- fT ^'»-( fT • (61) 

m— 1 m— 1 

The Hamiltonian -ff^ h commutes with j, J by construc- 
tion (the duality mapping preserves all algebraic rela- 
tions), meaning that no symmetry has been lost. Notice 
that j presents a highly non-local structure in terms of 
1 JJ II* 



Further Results on the Algebra of a Single Quantum 
Rigid Rotator 

Next it is shown that the operators j a defined in Eq. 
(42) satisfy the relations listed in Eqs. (45) and (47). The 
first step is to introduce the adjoint representation of the 
SU (2) Lie algebra via its double-covering homorphism R 
to 50(3), U i-> R(U), defined implicitly by 

U.a a .U* = J2<rbR(U) b a . (62) 

b 

Since tr(a a .a b )/2 = S ab , R{U) reads 

R{U) b a =tr(U.a a .U*.a b )/2. (63) 

It follows that Eq. (42) can be rewritten as j a = 
£ b i?(C/*)\j b . 



From Eq. (63), 

[J a ,R(U*) b c ] = ie abd R(U*) d c . (64) 

Also, \j a ,U] = ±J2 b U.a b R(U*) b a = \a a .U, where the 
last equality follows from Eq. (62) ( the conjugate re- 
lation [j a ,U*] = —U*.a a /2 follows in the same way). 
Combining this last result with Eq. (64) gives [j b , J a ] = 
iY. c ,d(- € adc + e cad )R{U*) c b J d = 0, and 

3'3 = EGC J bR{U*) b „)(£ R(U*) C a J c ) (65) 

a b c 

+ EtE^*)" - J *])(E R ( u *) c M = jj, 

a b c 

where the homomorphism property of R was used to sim- 

To check the commutator \j a ,jb], direct computation 
gives 

\j a> jb] = -te dec R(U*) d a R(U*) e b J c . (66) 

Since R(U) € 5*0(3), det R(U) = 1. It follows that 
e dec R{U*) d a R{U*) e b = e abf R{U*) c f . Then Eq. (66) 
simplifies to read 

[ja,jb] = -itabf R(U*) C f J c = -ie abc j c . (67) 

c 

It is only left to show that j\ = 
\Y. b {[J b MU*.<J a .U.<j b )\ + tr(U*.a a .U.a b )J b ) = j a . 
The commutator vanishes by virtue of Eq. (64) . 



* Electronic address: ecobaner@indiana.edu 
[1] H. Nishimori and G. Ortiz, Elements of Phase Transitions 

and Critical Phenomena (Oxford University Press, 2010). 
[2] E. Cobanera, G. Ortiz, and Z. Nussinov, Adv. in Phys., 

60, 679 (2011). 

[3] E. Fradkin and L. Susskind, Phys. Rev. D, 17, 2637 (1978). 

[4] J. B. Kogut, Rev. Mod. Phys., 51, 659 (1979). 

[5] J. K. L. Susskind, Phys. Rev. D, 11, 395 (1975). 

[6] J. Rosenberg, in Operator Algebras, Quantization and 
Noncommutative Geometry, A Centennial Celebration 
Honoring John von Neumann and Marshall H. Stone, 
Contemporary Mathematics, Vol. 365 (American Mathe- 
matical Society, http://www.ams.org/, 2004) pp. 331-354. 



